Experimental evidence for directed percolation in spatiotemporal intermittency 
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A new experimental system showing a transition to spatiotemporal intermittency is presented. 
It consists of a ring of hundred oscillating ferrofluidic spikes. The measured critical exponents of 
the system agree with those obtained from a theoretical model of directed percolation. 

PACS: 05.45-a, 47.54+r, 47.52 



Spatiotemporal Intermittency (STI) is characterized 
by patches of ordered and disordered states fluctuating 
stochastically in space and time. It is often a precursor of 
chaos in spatially extended systems In 1986 Pomeau 
suggested that the onset of chaos via STI might be 
analogous to directed percolation (DP) Q. Such pro- 
cesses are modeled as a probabilistic cellular automaton 
with two states per site, associated with the laminar and 



chaotic patches in the case of STI, and predict some uni- 
versal properties of STI. In particular, the fraction of 
chaotic domains is expected to grow with a power law 
e^, where e measures the distance from threshold. The 
mean laminar length decreases with e~^% as does the 
mean laminar time with e~^% and the critical distribu- 
tion is determined by H. 



TABLE I. Experiments and results. 



Authors 


Year 


Experiment 


Size 


T,{s) 


Geometry 








Ms 


Ciliberto et al. ffe 


1 1988 


RB-convection 


20 


10 


annular 




0.5 




1.9 ±0.1 


Daviaud et al. M 


1 1990 


RB-convection 


40 


2 


linear 


0.3 ±0.05 


0.5 ±0.05 


0.5 ±0.05 


1.6 ±0.2 








30 


2 


annular 




0.5 


0.5 


1.7 ±0.1 


Michalland et al. 


[| 1993 


viscous fingering 


40 


1.5 


linear 


0.45 ± 0.05 


0.5 




0.63 ± 0.02 


Willaime et al. | 


I 1993 


line of vortices 


15 


5 


linear 






0.5 




Degen et al. |10|l 


1996 


Taylor- Dean 


20 (90) 


1.5 


linear 


1.30 ±0.26 


^ 0.64 


^ 0.73 


1.67 ±0.14 


Colovas et al.~Ml 


1 1997 


Taylor- Couette 


30 (70) 


0.5 


linear 










Bottin et al. 


1997 


plane Couette 






linear 










Vallette et al. |13 


1 1997 


fluid fronts 


40 


0.5 


linear 










Jensen (theory) [ 


1 1999 


directed percolation 








0.276486(8) 


1.096854(4) 


1.733847(6) 


1.748 


present paper 




ferrofluidic spikes 


108 


0.08 


annular 


0.3 ±0.05 


1.2 ±0.1 


0.7 ±0.05 


1.7 ±0.05 
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FIG. 1. Scetch of the experimental setup. 



Some experimental tests of the conjecture in quasi 
one-dimensional systems have been made - ||l^. A 
short summary of these and the relevant exponents for 
the comparison with the DP-model are given in Table |. 
However the statement "...there is still no experiment 
where the critical behaviour of DP was seen" by Grass- 
berger still seems to be true. Thus in this paper a 
new experimental approach to this old problem is pre- 
sented. Motivated by the fact that a single peak of fer- 
rofluid shows chaotic oscillations under external driving 
of a magnetic field [|l^,|l^, we introduce here a system 
where about 100 of these oscillating peaks are coupled 
by magnetic and hydrodynamic interaction. They ex- 
hibit changes in peak height of about 10% and variations 
in wavelength A of about 50%. This system is advanta- 
geous because of its short response times and the easy 
control of the excitation. 

The setup of the experiment consists of a cylindrical 
electromagnet with a sharp edge (Fig. El). The magnetic 
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fluid is trapped by the inhomogenous magnetic field at 
this edge of the magnet icahy soft iron core. In that way 
the 40 mm diameter of the pole shoe supports a ring of 
up to 110 spikes of magnetic fluid as indicated by the pic- 
ture in Fig.|l|. The ring of spikes is recorded with a CCD- 
camera mounted above the pole shoe. The electromagnet 
contains a bias coil and an excitation coil. The bias coil 
is provided with a direct current of / = 1.0 A to keep the 
magnetic fluid in its place. This excitation coil is driven 
by an alternating current, phase- locked with the camera 
frequency, providing a stroboscopic jitter free recording 
on long timescales. The fluid is temperature controlled to 
12.5 ±0.03 °C by cooling the pole shoe, and encapsulated 
in a glass container to provide longterm stability. 




position X 

FIG. 2. Space time plots of the different states of the sys- 
tem, (a) lex = 2.8 A: laminar state; (b) lex = 2.8 A: spa- 
tiotemporal intermittency; (c) lex = 3.6 A: chaotic state. 500 
excitation periods are shown. The position is normalized over 
the size of the ring. 

The spatiotemporal behaviour is investigated by ob- 
serving the ring with the CCD-camera. To extract the 
wavelength and amplitude of the spikes as a function of 
space, we deflne a ring of interest around the center of 
the pole shoe which is covering the ring of spikes. It is 
divided into 1024 segments. The average of the grey val- 



ues within each segment represents the amplitude. This 
reduction of the 2D image to a single line scan can be 
done in real time with a frequency of 12.5 Hz. 




position X 

FIG. 3. x-t-plots at different states of the data processing. 

(a) raw data; (b) same section as in (a) after transforming the 
data to complex numbers. White corresponds to large local 
wavelengths, black to small ones, (c) final step: the relative 
change of local wavelength over two subsequent periods 

is calculated and binarized. The black areas are defined as 
regular (^A < 0.01), the white are chaotic (^A > 0.01). 

In Fig.^ 500 of such scans are shown in a space 
time plot, where dark regions correspond to high am- 
plitudes. The driving frequency of the excitation fleld is 
fex = 12.5 Hz, which is the same as the recording fre- 
quency. The period r = l/fex is used to scale the time. 
Due to this stroboscopic recording the regular oscilla- 
tions of the spikes cannot be seen. At low excitation 
amplitudes lex <C 3.0 A the system is completely regular 
(Fig.|| (a)) showing 108 spikes. Slight spatial variations 
of the wavelength A remain constant in time. In Fig.|| 

(b) at lex = 3.0 A irregular fluctuations are apparent. 

A fairly clear distinction between regular and irregular 
domains can be made in this image even by naked eye, 
which we consider as a manifestation of STL Further in- 
crease of lex leads to spreading of the irregular domains 
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engulfing the regular regions until finally the whole sys- 
tem is chaotic (Fig. ^ (c)). 

For a clearer distinction between regular and irregu- 
lar domains (Fig.|| (a)) we measure the local wavelength 
A(x,t) J|l7[, which is displayed using a grey scale plot 
in Fig.H(b). For the binarization the relative change 
SX = |At+i — At|/A^ of the local wavelength is used: 
Changes in ^A which are larger than 0.01 are counted as 
irregular, smaller changes belong to regular domains. In 
Fig.|| (c) the binarized values of SX are indicated. The 
quantitave analysis of STI described below is based on 
these binarized information. 

As an order parameter for STI we take the time- 
averaged chaotic fraction 7, which is the ratio of chaotic 
regions to the whole area of the system. Its variation 
with the control parameter I^x is shown in Fig.|^. Here 
the excitation amplitude is switched from I ex = 4.05 A to 
the value displayed on the x-axis. 7 is determined after a 
waiting time of ~ 1800r, when the transients have died 
out. Each data point represents a mean value of six in- 
dependent runs, which include six refills of the apparatus 
with fresh fluid. 
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excitation amplitude 7,^ (A) 
FIG. 4. The mean chaotic fraction 7 vs. excitation am- 
plitude lex- The solid line is a power law fit. The errorbars 
represent the statistical errors of the six measurements. 



Close to the onset of STI the mean chaotic fraction 
is expected to grow with a power law 7 ~ {lex — ^cY - 
The solid line in Fig.|4| is a fit to our data, using /c, 
and an offset representing background noise as ad- 
justable parameters. The threshold value determined 
in this way is /c = 3.0 ± 0.03 A and the exponent 
(3 = 0.3 ± 0.05 is in agreement with the theoretical ex- 
pectation (3 = 0.276486(8) g. 




control parameter e 
FIG. 5. (a) The mean chaotic fraction 7 vs. control pa- 
rameter e. (b) The mean laminar expansion in space i and 
in time T vs. control parameter e. The solid lines are power 
law fits. 



An alternative representation of the same data is pre- 
sented in Fig.^ (a) in a log-log plot. Here the normalized 
control parameter e = lex/Ic — ^ serves to define the 
x-axis. In Fig. I (b) the mean laminar length which 
is the average of the spatial expansion of all regular do- 
mains at a given e, is presented. It is expected to grow 
with a power law i ~ e~^^ and the fit yields an expo- 
nent Us = 1.2 ± 0.1, which is represented by the solid 
line. Only data in the range of 0.03 < e < 0.1 are taken 
into account. For smaller e a finite size effect is obvious: 
i must be smaller than the system size 1. For larger e 
the system is no longer intermittent, but rather chaotic. 
The same argument holds for the mean laminar time T, 
which is the average of the life time of all regular do- 
mains (Fig.^(b)). A similar power law behaviour can be 
seen with an exponent Ut = 0.7 ± 0.05. The exponent 
Us is in good comparison to theoretical expectations of 
Us = 1.096854(4), while Ut differs from the theoretical 
value ut = 1.733847(6). 

In Fig. H the distribution of the laminar domain length 
/ for e = is presented. The solid line represents a fit 
p{l) ~ with /is = 1.7 ± 0.05, in agreement with the 
theoretical value jUs = 1.734. Within the experimental 
errors, the theoretically predicted relation 

Ma = 2 - ^ (1) 
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is thus fulfilled. 
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FIG. 6. Distribution of the laminar domain length / for 
e = 0. The solid line is a power law fit. To suppress the sta- 
tistical fluctuations the distribution is logarithmical binned. 

The distribution for e > as shown in Fig. ^ cannot 
be described by a simple power law. A more complicated 
distribution function has been suggested 

p(0 = (Ar^ + 5)e-^/^^— ^ (2) 

in Ref. . The solid line is a fit to this empirical distri- 
bution function for e = 0.019 with Idecay = 0.17. It shows 
clearly that the power law is now replaced by a function 
more reminiscent of an exponential decay. 
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laminar length I 

FIG. 7. Distribution of the laminar domain length / for 
e = 0.019. The solid line is from eq. 

In summary, we have measured the critical exponents 
/3, z^s, vt and jig for the mean chaotic fraction 7, the 
mean laminar length ^, the mean laminar time T and the 
length distribution function. Three of these four expo- 
nents agree with the theoretical expectation derived from 
a DP model and thus fulfill the relation (|l|). When con- 
sidering the fact that the theory is applicable only near 
Ic and that our apparatus is influenced by finite size ef- 
fects, the results support the idea that our system shows 



a DP like transition to STI. Further investigations will 
thus focus on the influence of the system size. 
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